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1 /INTRODUCTION 


r 


This monograph discusses the lateral model development and modal calculations 
for the cylindrical space vehicle system and also systems employing clustered 
tanks. Relative importance of physical characteristics will be discussed as well 
as methods used by the industry for the solution of modal parameters. Primary 
attention is focused on parameters important in control and stability analyses for 
which the system frequency of interest is generally below twenty cycles per second 
and quite often below ten cycles per second. Application for loads analysis follows 
the same principals outlined herein but may require more detailed representation 
in areas where loads or deflections can be critical. The models described will be 
concerned with gross vehicle lateral motions and will not include the shell radial 
deflection modes which can be significant for local effects analysis but generally 
do not contribute to gross vehicle control and stability. Related monographs describe 
torsional, longitudinal, and sloshing models. Stability and loads analytical methods 
using these models and modal parameters are also subjects of other monographs. 

The general approach for dynamic solutions involving large systems is to develop 
a mathematical model describing the system's mass and structure, calculate its normal 
modes of vibration, and then, using normal mode theory, apply the external forces 
and couple in the control system to obtain total response. The dynamic analysis 
is then only as accurate as provided by the mathematical model representing the 
space vehicle system; therefore, development of these models is of major importance 
in dynamic analysis. Also, since these models are idealizations and approximations 
of the real system, the experience of an analyst in deciding which elements are 
dominant contributes greatly to the successful representation of the system. 

Representation of the space vehicle system in the lateral direction is accomplished 
by a series of lumped masses connected by elastic beams. By successive refinements, 
such as branch beams to include multiple load paths, concentrated masses attached 
to the beam by translational and rotational springs, or other possible independent 
structures or components, the model can simulate all significant motions. The 
mathematical model and numerical techniques to obtain lateral normal modes for 
space vehicle systems comprised of tandem cylindrical structures have been 
extensively developed by aerospace industries. 

The addition of peripheral tanks to a center core increases the complexity by 
three mechanisms; (1) the attachments between tanks require additional boundary 
conditions: (2) the attaching structure is generally complex, and (3) the peripheral 
tanks generally destroy axis-symmetry and results in multi-directional normal modes. 
The present techniques employed with these clustered boosters are essentially the 
same as those of the simple cylindrical structure; however, the complexity dictates 
use of more points to describe the system and the analysis can easily exceed computer 
storage capacity. This leads to a compromise between accuracy of results and 
efficiency and capability of present computer facilities. 
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2 /STATE -OF -THE -ART 


The original work in the study of the phenomenon of structural dynamic behavior 
was primarily restricted to simple spring-mass systems and severly limited in 
its coverage of continuous systems to uniform beams and plates. As structural 
complexity increased, consideration of practicality rendered these solutions 
unjustified. The analysis of larger systems required vast amounts of time and 
effort and, since all calculations were done by hand, numerical accuracy became 
a problem of major concern. The most practical solution to a structural dynamics 
problem generally turned out to be a healthy factor of safety applied to structural 
capacity. 

With the advent of aircraft, a large safety factor could no longer provide the 
means of accomodating dynamic effects. The development of machine calculators 
enabled the analyst to develop more complex dynamic models and associated 
solution techniques to provide answers to the problems of stability and flutter in 
aircraft design. Still further refinement in modeling techniques has been permitted 
by the introduction of electronic computing systems. The capacity of these 
computing systems and rapid computation cycles has resulted in the development 
of more numerical techniques utilizing smaller iteration intervals with a 
resultant increases in solution accuracy. 

Approximate methods for solving the equations resultant from model simulations 
have been developed; a summary (from Reference 1) in tabular form of the 
representative accuracy of some of these techniques is presented in Table 1. 

To illustrate the ability of lumped parameter systems to represent continuous 
systems, consider a uniform simply supported beam. The exact expressions for 
the natural frequencies and their associated mode shapes are well known: 

«•-<&)* f 

ip n = Sin np 

where n is mode number, 1 is the length of the beam, and H is the mass per 
unit length. Assuming the following physical data: 

1 = 120 inches (3.05 meters) 

9 2 13 2 

El = 8 x 10 lb-in (2.33 x 10 gr cm ) 

M = 400 lb /inch (7. 14 x 10^ gr/cm) 
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the resultant first three natural frequencies are: 
«! = 3626.5; 0^= 60,22; f=9.58cps 
W® = 58024.; U) s = 240. 88; f = 38. 34 cps 

60^ = 293748. ; U> 3 = 541.98; f = 86.26 cps 


Table 1. Comparison of Accuracy of Various 
Methods of Approximate Solution 


Method 

Number of 
Masses 

CO 

1 

% Error 

W 2 

% Error 

Rayleigh-Ritz 
2 -Polynomial 

10 

2.5372 

2.38 

12.9085 

42.0 

Rayleigh-Ritz 

3-Polynomial 

10 

2.4894 

0.45 

9.9379 

9.34 

Modified 

Rayleigh-Ritz 

2-Polynomial 

10 

2.4769 

-0.05 

9.4832 

4.34 

Collo.cation 

11 

2.4760 

-0.09 

8.9207 

-1.85 

2 Station 
Functions 

11 

2.4773 

-0.04 

9. 0012 

-0.96 

Iterated 

Rayleigh-Ritz 

Polynomial 

(Galerkin) 

10 

2.4769 

-0.05 

9.4832 

4.34 

Holzer- 
My kies tad 

10 

2.4776 

-0.02 

9.1513 

0.69 

Direct 

Iteration 

5 

2.4761 

-0.09 

9.3457 

2.82 
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The associated mode shapes $ 1 , $2 , and 0a are sine waves of one, two and three 
halfwaves, respectively. 

Simulating this same simply supported beam by a set of equal lumped masses 
spaced equidistant along the beam establishes a discrete system which may be 
solved in any number of different ways. Trial cases have been analyzed using 
10, 20, and 30 lumped masses. Comparison of the resultant natural frequencies 
can be seen in Table 2. 

The accuracy of the solution of the discrete system varies with the number of 
degrees of freedom given to it; note, however, that the relationship is not direct; 
the accuracy requirements of the task influence the decision as to how many 
degrees of freedom can be justifiably assigned to the model. 


Table 2. Frequency Comparison 


Exact Solution and Discrete System Solution for Simply 
Supported Uniform Beam 


^\§ystem 

Number 

Continuous 

(Exact 

Solution) 

10 -mass 

lumped 

parameter 

20-mass 

lumped 

parameter 

30-mass 

lumped 

parameter 

1 

9.58 

9.58 

9.58 

9. 58 

2 

38.34 

38.33 

38.34 

38.34 

3 

86.26 

86,20 

86.26 

86.26 


(Frequencies are in cycles per second, ) 


Application of the lumped parameter techniques to space vehicle systems has 
been used extensively by the aerospace industry. Invariably, the propellant 
tanks and all cylindrical structures are described as beams with propellant and 
structural mass distributed. Local stiffnesses are represented as equivalent 
translational and rotational springs or by experimentally or analytically 
determined influence coefficient matrices. 

There have been several attempts to verify the dynamic characteristics 
of a launch vehicle as indicated by the analytical model with those obtained 
by test on the prototype itself, A comparison of experimental and analytical 
modes for an Atlas with a long, slender payload is given in Reference 2. 

These results show that frequencies for the first four modes are within 10% 
accuracy. The mode shapes are fairly well represented in the first four 
bending modes although displacement amplitudes were found to inorease in 
error with mode number, Using the generalized mass as a measure of mode 
displacement error these are 4%, 28.2%, 177%, and 20% for the first 4 modes. 
The fifth mode shape was not well represented although frequency comparison 
was within 15%, 
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A similar comparison was made for the Minuteman solid ICBM in Reference 
3. These results give a 13% difference for frequency in the first four modes 
and slightly larger errors in mode shape. After a correction for adapter joint 
flexibility the maximum error in frequency of the first four modes was 6% 
with corresponding improvement in mode shape. 

Vibration test results of an Atlas-Agena-OAO vehicle, which is a more 
complex structure than the previously mentioned vehicles, are presented in 
Reference 4. Here, the first mode compares reasonably well with the 
analytical modes; however, the higher modes were not representative of the 
analytical modes. Detailed examination of local structure revealed possible 
amplitude dependent component stiffnesses. When this characteristic was 
taken into consideration, the comparison between analytical and test mode 
shapes was improved to a satisfactory level, 

Several important factors can be obtained from the results of these 
comparisons: (1) the idealization of the tank as a one -dimensional beam 
appears to be adequate ; {2) much attention should be focused on local 
structure such as connections; (3) the numerical procedures presently 
available are capable of solving for the primary modes. These studies 
also verify what might be inferred from the mathematical characteristics: 
just as the eigenvalues (frequencies) can be determined with greater accuracy 
than the eigenvectors (mode shapes) , so are the frequencies of the physical 
system more accurately predicted analytically than the corresponding mode 
shapes. Furthermore, with the numerical methods used to analyze large 
systems, the accuracy of the calculated characteristics decreases with 
increasing order of mode number, The redeeming feature is the fact 
that for most vehicles the overwhelming fraction of response is obtained 
from the first and sometimes second mode, Consequently, although a 
sizeable error may be Introduced in the higher modes, the net error in 
total vehicle response is small. As a matter of fact, in the oase of cylindrical 
vehicles, attempts to verify the mathematical representation in flight has 
succeeded in only identifying the first mode and second mode frequencies 
from the flight data, 

Analysis of clustered boosters is still in the formative stage, although 
much improved over the earliest attempts at analysis. The first efforts 
(Reference 5) attempted single beam models, but when tests on the vehicle 
were performed, the total inadequacies of the single beam analogy were made 
apparent. It was evident that multidimensional motion and individual tank 
motions must be defined analytically. The analysis procedure was revised 
to provide for such motions. The first attempts at such analysis produced 
results which were still significantly atvarianoe with test data. The primary 
cause was found in the evaluation of the stiffness of the complex structure. 

As structural definition Improved, so did the modal data. Mathematical 
difficulties imposed by excessive numbers of coordinates have been overoome 
by use of the modal synthesis technique (Seotlon 4. 3). Modal data thus 
obtained compares well with both static and flight test data. 
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In the case of Saturn, the frequency comparison between analysis and test is quite 
good, (Reference 6). Again the mode shapes do not agree as well, especially 
in the motion of peripheral tanks. The effect of these modal discrepancies on 
control and loads is not important. Techniques have been developed whereby 
these discrepancies are accounted for and reworked as they appear. These 
discrepancies also vary from one vehicle to another, which sometimes necessitates 
additional work to obtain satisfactory agreement between model and test. 

A variance with cylindrical tank vehicle experience is the composition of flight 
test response data of clustered tank vehicles . Such data from gyro and accelerometer 
traces indicates measurable contributions to vehicle response from several bending 
modes and higher complex modes. (See Reference 7). 
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3 /MODEL REQUIREMENTS AND RECOMMENDED PROCEDURES 

The solution of dynamics problems is approached by developing a set of governing 
equations through consideration of the condition of dynamic equilibrium or of 
energy relationships of the system under applied external forces. For an un- 
damped system the general form of these equations, in matrix notation, is 

r.M 1 (z } + [K 1 [z } = {F ] ( 1 ) 

where Cm Jis a matrix of masses, [K ] is a stiffness matrix, { z } and { % 3 
are the displacement and acceleration vectors, respectively, and [F } is 
a vector of external forces. 

A very useful characteristic of elastic systems is that they will respond or 
vibrate in natural orthogonal modes of displacement $n with circular 
frequencies co n . The total displacement of a system can then be expressed 
as a summation of individual natural mode displacements. 

This is given by 

U } = [*] [q] (2) 

where [ $ Jis a matrix of mode shapes and q n is the time dependent amplitude 
of mode n. Substituting equation (2) into equation (1) yields 

[M] [$> J U’ J + [K] [4-J [q } = { F}. (3) 


Premultiplying both sides of (3) by [.$ ] gives 

/ / 

[*1 [M][$]{d'} + [*] CkI [*] {q} = [40 [F}. <4) 

the requirements of orthogonality and harmonic motion of the natural modes 
provide the relationships 

[$> V [M ] [ <f> ] = fmj (5) 

and 

[*] # [K] [* J = tVJ CmJ . (6) 

Substituting (5) and (6) into (4), we obtain 

U } + r CO 3 ] C q } = tmf 1 [4 3 [F } -[m] -1 {Q}. (7) 
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Equation (7) is a set of n uncoupled equations in terms of qn, to n , the generalized 
mass, m n and the generalized force Q n , The solutions of these equations identify 
the time dependent values of qn which are then used in equation (2) to give complete 
system response. Detailed discussions, derivations, and proofs of the equations 
of motion, orthogonality of natural modes, and normal mode theory are given in 
References 8 to 11. 

The use of normal mode theory requires determination of these natural modes of 
vibration. If harmonic motion is assumed and the applied forces are equal to zero, 
then equation (1) can be written as 

- to 2 [M ] {z } + [K ] { z } = 0 


or 

Cz } = to 2 [K]" 1 [Ml fa}. 


( 8 ) 


Each of these equations is in a form suitable for solutions to obtain the orthogonal 
modes and their natural frequencies. Many numerical techniques have been 
developed to obtain these characteristics and several are discussed in Section 4.2. 

The above discussion presents the fundamental approach to structural dynamic 
response analyses. The mathematical model used to simulate the physical system 
is the basis for determining the natural modes and the solution of the dynamic 
problem. It is significant, therefore, to examine in detail the factors influencing 
its development. 

3. 1 LATERAL REPRESENTATION OF CYLINDRICAL LIQUID PROPELLANT 

VEHICLE. In most instances, the lateral dynamic characteristics of liquid 
rocket propelled space vehicle systems can be considered to be adequately repre- 
sented by simple one-dimensional beam theory. It is common practice, and 
certainly more convenient, to replace the continuous structure by a lumped parameter 
Idealization, In such an idealization, the analyst concentrates attention on those 
aspects of the system which are felt to be dominant (major masses, major structural 
elements, propellants). The disorete model is formed by oonoentrating mass at 
selected points along the beam, These points are ideally the o.g. of the distributed 
mass which is to be considered concentrated at that point. 

Elastic properties are expressed in lumped fashion as a set of flexibility co- 
efficients, Cy , or stiffness influence coefficients, Ku . These ooeffiolents have a 
physical significance in that Cm can be considered as the deflection of point 1 due to 
a unit load at j; and Kn equated to the force produced at point 1 due to a unit deflection 
at point j, if all coordinates other than j are temporarily restrained. (Flexibility and 


stiffness influence coeff cients are covered in more detail in Section 4. 1). 

The mathematical description of this discrete model is a set of simultaneous 
linear ordinary differential equations. Such equations lend themselves readily 
to matrix techniques and the utilization of digital computers permits the solutions 
of sets of equations too large for practical hand computation. 

The one-dimensional beam representation is the simplest lateral model and 
may not fulfill all necessary requirements for a specific problem. This can 
necessitate recognition of non- structural modes (sloshing), local response 
characteristics (engines) or multiple load paths not accounted for in the simple 
beam analogy. A further refinement of the model is then necessary. 

3.1.1 MASS AND ROTARY INERTIA. The distributed mass and inertia data 
must be lumped into discrete, point masses, the number of which determine the 
degrees of freedom given to the model . The number of mass stations is influenced 
by the number of bending modes to be calculated. 

It has been found that for one-dimensional beam bending models the required 
number of mass stations should be approximately ten times the number corresponding 
to the highest elastic bending mode to be calculated . For example, if three elastic 
bending modes are to be calculated, then approximately thirty mass stations are 
required to represent adequately the bending dynamics of the third mode. This 
criterion has been established empirically by calculating mode shape, frequency, and 
generalized mass corresponding to the first three elastic bending modes for typical 
vehicles configurations in which the numbers of mass stations used were successively 
increased from eighteen to forty. As expected, the acouraoy Increased as additional 
stations were utilized. However, it was observed that no further significant inorease 
in accuracy was achieved by using more than thirty mass stations. 

For a more complex model (such as one with branched beams) the above general 
rule may not be strictly applicable. For a branohed system, the general rule may 
be applied to the primary beam of the system, and masses lumped on the secondary 
branches in about the same distribution. It must be emphasized that as the model 
diverges from the single beam oonoept, performed mass lumping rules become less 
applicable and more reliance must be plaoed upon the experience of the analyist. 

Note that only rigid masses are to be inoluded in this distribution, that is, only 
those masses which can be considered to aot as an integral part of the unrestrained 
beam during its vibrations. It oannot be over-emphasized that items such as pumps, 
equipment pods, eto. , which are aotually, or simulated as, mounted elastically to the 
main structure may significantly alter the bending oharaoteristios of the higher 
frequency modes. 
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Whether or not such masses are to be treated as integral to the beam or as separate", 
elastically attached masses, depends upon: (1) whether or not the frequencies 
of the body modes to be computed are less than or greater than the mount frequencies 
of the discrete masses, and (2) whether or not these masses are great enough to 
materially affect the result. 

Accurate representation of the distributed mass at discrete points would require 
inclusion of the mass moment of inertia of the distributed mass at that point. 

With liquid propellant, the effective moment of inertia is not easily determined. 
Fortunately, these rotary inertias have only a small influence on the modal quantities 
of the first two bending modes and can be neglected as shown in the work of Reference 12. 
12. This work is a comparison of modal parameters for a typical launch vehicle 
for calculations with and without rotary inertia terms. (Rotary inertia was 
calculated assuming the liquid behaved as a solid). The differences in the first 
two modes were less than 10% for all parameters and less than 5% for most para- 
meters, In the third and fourth modes the differences were usually less than 20% 
with a few extreme differences up to 100%. These calculations are at the extremes 
of effective rotary inertia and the correct solution is probably somewhere between. 

A better representation could be obtained by using more mass stations rather than 
including rotary inertia. 

Because of the small effect in the lower bending modes and the uncertainty 
of effective inertia, it has been common practice to neglect this term. Also, this 
allows for much more efficient computation since this eliminates half the coordinates 
in the solution of characteristic equations. 

3.1.2 SLOSHING PROPELLANTS. Space vehicle system propellants constitute a 
large percentage of the total system weight. Part of this propellant can be con- 
sidered as rigid or distributed mass on the Idealized beam while a smaller portion 
must be allowed to slosh In the lateral model. This sloshing mass becomes more 
important in later flight times, when the sloshing mass becomes a sizeable pro- 
portion of total vehicle weight. The complex problem of control system coupling 
with multi-tank slosh masses precludes any simple rule of thumb criteria for 
sloshing stability. Experiences with some vehicles indicate major instability 
periods occur when the sloshing masses approach 20% of total vehicle mass. 

Several methods have been developed for the description of propellant sloshing 
modes and frequencies, The general approach, as related to lateral models, is 
to derive the hydrodynamio equations into a form suitable for a meohanioal analogy. 

It can be shown for a oylindrioal tank that if: (1) the tank walls are rigid, (2) the 
fluid is incompressible and lrrotational, and (3) only small disturbances are 
admitted, then pendulum or spring mass analogies oan be devised whioh will re- 
produce the oharaoteristios of the fundamental mode of sloshing oscillation. 

Attaching the equivalent spring-mass to the lateral model la easily aooomplished 
and is therefore preferred over the pendulum. The seoond and higher sloshing 
modes are not generally considered beoause the magnitude of lateral foroe 
contribution from these modes decreases rapidly with increasing order; further- 
more, test experience indicates a great deal of turbulent mixing occurs and 
damping effeots are greater in higher modes. 
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Techniques for deriving the mechanical analogies for slosh and their 
limitations are the subject of another monograph. 

3.1.3 ENGINE REPRESENTATION. Thrust-vector control of liquid- 
propellant vehicles is generally maintained by gimballing the rocket engines. 
Since the entire engine is gimballed rather than just the thrust vector, this 
gimballing action will cause inertial forces as well as thrust forces to act 
on the missile body. These inertial forces are appreciable, and their 
lateral components will equal those of the thrust forces when the engine is 
gimballed sinusoidally at a particular frequency; at higher frequencies 

they can exceed thrust forces. The thrust vector displacement is determined 
by two effects (1) the displacement contained within the elastic mode 
(including the flexibility of the engine mounting and actuator structure) while 
the servo positioning system is locked and (2) the additional degree of 
rotational freedom added to represent the motion accompanying action of the 
positioning servo. For an adequate representation in an analysis including 
the engine and control systems it is necessary to include engine characteristics 
in the lateral model. The only exception arises when the model frequencies 
are well below those of the engine mounting system. 

The engine is incorporated into the lateral model by attaching a mass and 
moment of inertia at the appropriate location on the one -dimensional beam. 
Since the engine itself is quite rigid, the only elasticity normally considered is 
the mounting structure and actuator system. This structure is generally 
complex and test data is often required for proper simulation. One such 
test would be a vibration test to determine the resonant frequency of the engine 
on its mounts and using these results obtain the equivalent rotational spring 
connecting the engine to the vehicle. This primary frequency is often low 
enough to fall within the range of the lower vehicle bending frequencies and 
as a result could have a significant effect on bending stability. 

Other means of obtaining thrust vector control are used and require engine 
representation in the lateral model. Vehicles using a combination of fixed 
and gimballing or swivel engines require simulation of the fixed engines as 
well as the engines with freedom to represent their displacement in the 
vehicle elastic modes. Control concepts such as movable nozzles or stream 
deflection, involve little or no additional mass motion and, therefore, only 
the fixed engine representation would have a significant effect on the lateral 
modes. 

3.1.4 BRANCH BEAMS, Frequently, the vehicle construction will be such 
that major portions are cantilevered within another structure or are connected 
through different load paths. Examples are: payloads enveloped by protective 
fairings, engine compartments of upper stage vehicles suspended in the inter- 
stage adapter well, or multi-engine vehicles having independent load paths for 
each engine - such as a center engine supported on the tank cone and peripheral 
engines mounted to the cylindrical structure of the vehicle. 
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Such conditions are illustrated in Figure 1. Realistic representations of 
these arrangements are required, not only for true definition of gross vehicle 
response, but also to investigate possible interference between parts. These 
multiple paths can be accounted for by appropriate branch beams from the major 
planar beam. So long as the analysis is restrained to one dimensional motion, 
there is no significant added complexity introduced by the branch beams since 
the compatibility relationships at the junction 



points can be easily satisfied. Note the model in Figure 1, Branch beams can 
be attached in two ways: by secondary beam elements, as is done for the payload 
fairing and the upper stage engine structure ; or by concentrating elasticity in 
lateral linear and angular springs, as is done for the external engines of the 
booster. The only mathematical consideration involving in choosing which 
analogy is more appropriate is that the beam element will influence more 
elements of the flexibility (or stiffness) matrix than springs will, due to the 
off diagonal cross-coupling terms. 

Generally, these branoh beam oonditions are encountered with relatively 
small masses and do not alter the gross vehicle modes significantly except 
that an additional mode is accounted for where the two branches are out of 
phase but at nearly the same frequency as the in phase mode. When these branch 
beams involve engine displacements, they can be significant for control and 
stability analyses. For other portions of the vehicle they should be included 
to obtain proper load distribution and clearance envelopes. 
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3. 1. 5 LOCAL STRUCTURE EFFECTS. One of the major difficulties en- 
countered in describing a vehicle is the effect of local structures such as 
joints between the interstage adapters and vehicle stages, trusses on which 
payload or engines are mounted, or play in joints such as engine gimbal 
blocks when the engine is not under thrust. In the case of an adapter joint, 
its stiffness may vary under compression or tension. Although these joints 
are usually in compression, it Is possible during the period of maximum 
aerodynamic loading for the combination of axial and bending loads to cause 
one side of the cylinder to be in tension. Depending on the characteri sties of 
the joint, this could lead to a significant error in frequency and mode shape. 

The variation in stiffness of the joint under these conditions is difficult to 
determine accurately by analysis and usually test verification is necessary to 
determine the significance of this effect. Once these values are obtained, 
they can be substituted in the model and used for the modal calculations. 

Similar problems can exist for the local structure supporting engines 
or payloads since these structures are often redundant, carrying loads to 
a flexible shell. It is possible to obtain these influence coefficients analytically, 
but a final check with test results is advisable. The free play occasionally 
found in connections such as an upper stage engine gimbal (when not under 
thrust) during first stage flight is random and difficult to represent. These 
can produce some low frequency pendulum or inverted pendulum modes of 
significance if the mass involved is appreciable. To determine such effects 
some crude pendulum-spring mass analogies can be used to establish whether 
or not further consideration is necessary. 

The above are a few local effects to be examined in construction of the 
lateral model. In general, joints that carry significant loads or components 
of sizeable mass should be examined in some detail to establish the degree 
of representation required in the lateral model. 

3.1.6 LOCAL NON-LINEARITIES. If major non-linearities exist, the system 
and its response cannot be described correctly with conventional normal mode 
analysis techniques. The effect of a separation joint possessing non-linear 
bending stiffness was investigated through the use of quasi -normal modes and 
a Rayleigh - Rltz analysis in the work of Reference 12. In the analysis, the 
assumed mode shapes are those of the vehicle having an infinitely stiff 
separation joint plus one additional mode having a single concentrated non- 
linear rotational spring located at the separation point with the remainder 
of the vehicle considered as rigid. The LaGrange equations produced 
simultaneous equations in the normal mode coordinates with inertial 
coupling between the orthogonal elastic modes and the non-linear spring 
mode. 
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To the equations of motion developed with the above techniques were added 
the control sensors, engine representation, and control system representation 
for a bending stability analysis. Because of the non-linearities in both the 
vehicle structure and the engine actuators the solution was obtained with an 
analog computer. The study presents an approach for solving problems in 
structures with non-linearities using models modified to account for local 
peculiarities. 

3.1.7 TEMPERATURE. Space vehicle systems primary structure are 
subjected to change in temperature of hundreds of degrees varying from 
cryogenic temperature to the extreme elevated temperature resulting from 
aerodynamic heating. This increase in temperature causes a reduction in 
the material moduli which in turn leads to a small reduction in frequencies 
and altered mode shapes. Temperature considerations are unimportant 
until after the period of maximum aerodynamic pressure and then only for 
certain portions of the vehicle. Since the period of maximum heating 
occurs after the period of maximum disturbance and only affects parts of the 
structure, its significance is greatly reduced. The heating of various 
portions of the vehicle can be predicted within tolerances necessary for 
modal analyses to establish the resultant variation in modal parameters. 

It is a condition which should be examined if bending stability is marginal. 

3. 1. 8 AXIAL LOAD, Axial loads caused by longitudinal acceleration of 
several g's during flight will cause a slight decrease in bending mode 
frequency through two mechanisms: (1) the effect of axial load on beam vibration, 
and (2) the reduction in equivalent skin on stringer-skin structure. The 
first effect can be represented analytically and the second can be included 
after calculating or obtaining the equivalent skin from empirical data. 

The total effect of axial loads is generally very small and in nearly all 
cases can be ignored. In cases where the control system Is marginally 
stable this is one of the several minor effects that must be evaluated. 

The study in Reference 13 shows that the effect on frequency for 
representative booster structure is less than 1% when shear stresses are 
small. 

3.1.9 LATERAL-TORSIONAL-LONGITUDINAL COUPLING. The typical 
axis symmetric cylindrical space vehicle is analyzed as if lateral , torsional , 
and longitudinal motion are not ooupled. Actually, these vehicles are not 
completely symmetric and a possible ooupllng mechanism, however slight, 
can always be found. The importance of this coupling oan vary greatly 
from vehiole to vehicle and even if it is known to exist from flight or 
experimental data, the coupling mechanism is difficult to identify. These 
coupling problems often ocour when the modal frequencies of two modes, 
say, one lateral and one torsional, are very olose together. Then a very 
small coupling mechanism, suoh as c.g. offset from the supposed line of 
symmetry, can result in coupled motion. 
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A comparions of the frequencies of the modes in the three directions 
should be made to determine the existence of modes of nearly equal frequency. 

If such a condition exists , it is necessary to examine the condition under 
which this may cause a significant problem. As an example, if a limit cycle 
can occur due to sloshing, could this cause excitation of a critical torsional 
mode at this same frequency or sub-harmonic ? In most instances of coupling 
of this type, a periodic forcing function is necessary to transmit the energy 
from one direction to another. 

Cylindrical vehicles with unsymmetric upper stages or payloads of large 
mass can cause coupling in the various directions in the low frequency modes. 
The model (and analysis) then become complicated and approaches that of the 
clustered boosters. Representation of this configuration requires detailed 
description in the unsymmetric stages and proceeding with analysis as 
described later for clustered boosters. Preliminary work would indicate 
the degree of sophistication to be used for adequate representation for 
stability and loads analysis, 

3.1.10 DAMPING EFFECTS. Structural dissipative (damping) forces exist 
in the vibrating structure as a result of material strain hysteresis and coulomb 
friction in structural joints. The nature of these damping effects is obscure 
and does not lend itself to analysis other than an approximate empirical 
treatment, by which the gross effect of these scattered dissipative mechanisms 
is represented as an equivalent viscous damping, added to each mode as 
appropriate. The damping is thus assumed to produce no coupling between 
modes. While this mechanization is not entirely realistic, it is justified 
by the following observations: 

a. The actual damping is very low and is found by test to produce little 
coupling. Thus, nearly pure normal modes of a system may be 
excited and the system observed to decay almost harmonically. 

The indication given is that velocity dependent coupling is very small. 

b. If an attempt is made to show a velocity dependent coupling, the 
coefficient would have to be determined experimentally. Since 
the direct damping coefficient is itself difficult enough to measure 
it is clear that the accuracy of a study can not be increased by the 
introduction of still more suspect data. 

The structural damping force is a function of the deflection of the 
generalized coordinate of the mode but in phase with the velocity of the 
generalized coordinate of that mode. To treat this damping as a viscous 
damping requires that the mode oscillate in a quasi -harmonic manner. 

This damping force may then be expressed as a damping factor, £ n> 
where 2£ n o) n 4n is tee internal damping force of the nth mode per unit 
generalized mass. 
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Fluid propellant damping forces result from the dissipative nature of a 
viscous fluid undergoing shear. Although there are some approximate 
methods for calculating damping forces, these forces are most commonly 
arrived at by experimental testing of the actual tank, in the case of small 
missiles, and a model tank In the case of large missiles. These forces may 
be represented as a propellant damping factor, ^ , in the expression 

2 ^l w l4l which is the damping force per unit sloshing mass ahd cii is 
the lateral velocity of the 1th sloshing mass. 

Equation (7), with damping included, becomes 

tfn + 2 in “n q n + w n % = m n X $n (9) 

Methods for obtaining values of £ n from vibration test date are given 
in the monograph covering that subject. The dissipative forces associated 
with sloshing are covered in the sloshing model monograph. 

3.2 ADDING COMPONENTS USING MODE SYNTHESIS. Frequently it is 
desirable to make a parameter study to determine the effect on vehicle 
response resulting from changes in the characteristics of a specific area 
or component, e.g. a sloshing mass or engine system. Rather than make 
several analyses of the system changing but a fraction of the parameters 
each time, the vibration characteristics of the system excluding the specific 
varying parameter may be calculated, and then modified by coupling the 
parameter back in through the mode synthesis technique (discussed in 
Section 4,3), 

For reasons developed more fully in Section 4.3, the mode synthesis 
approach may result in a loss of accuracy. The analysis that considers 
the most information about the system will be the most accurate, The 
use of many modes in the mode synthesis technique will give theoretically 
more accurate results than using a minimum number of modes. This 
aspect is one which must be handled by discretion born of experience. As 
an example, in calculating the bending modes with sloshing propellants, 
three alternatives are available : 

1. Include the sloshing mass as an attached spring-mass to the beam 
in the modes calculation. 

2, Assume the sloshing mass is included in the "rigid" propellants 
for modes calculation. The sloshing is then included through 
mode synthesis by adding the single spring mass mode and 
subtracting the sloshing mass effects from the bending modes. 

This requires both inertial and elastic coupling in the synthesis. 
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3. Assume the sloshing mass can be eliminated in the "rigid" pro- 
pellants for modes calculations. The sloshing is then included 
through mode synthesis by adding the simple spring mass mode. 

This requires only elastic coupling in the synthesis. 

These alternatives are listed in order of accuracy of end results 
following the general rule stated previously. There are many examples 
where there would be little if any degradation of accuracy. As an 
example, consider a vehicle with first bending frequency of 5 cps 
and first slosh frequency of 1 cps, then the sloshing is essentially 
uncoupled from the elastic modes. If the slosh frequency were 4 
cps, then considerable coupling is possible. 

Although sloshing was used as an example, the same is true for 
any representation of this type, i.e. , engines, payloads, sloshing 
propellants, or any other large component, or specific parameter under 
investigation. 

3.3 CORRECTING MODEL BASED ON TEST RESULTS. The final 
verification of analytical techniques is a comparison with experi- 
mental data. Perfect comparisons are indeed exceptions, since 
both the analytical model and experimental model are approximations 
to some extent. The analytical approximations have been discussed. 

The major experimental approximations are centered around suspension 
system effects and vehicle modifications required to accommodate the 
suspension system. No general rule can be made to obtain better agreement 
between test and analysis. Careful examination of the data and the 
structure will probably indicate several areas where the representation 
is inadequate or does not define the test specimen. Possible causes of 
differences are: 

(1) Effects of suspension system on test environment. 

(2) Stiffness of joints or trusses. 

(3) Assumed planes of symmetry incorrect. 

(4) Effeot of large component* , suoh a* engine*. 

(5) Experimental mode* may be Impure, l.e, , not orthogonal or 
inolude parts of other mode*. 

(6) Effeot* of rotary inertia. 

The work of Referenoe 14 presents a method for obtaining the flexibility 
matrix from experimental mode data. The prooedure orthogonal ize* the 
experimental modes, using an analytical mass distribution, and then derives 
the flexibility matrix of the struoture. This method oan be useful if complete 
and aoourate experimental data Is obtained for a system diffloult to model. 

It oan also be used to looate possible dlsorepanoies between analytioal and 
experimental results. 
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3. 4 SOLID BOOSTERS. The solid propellant grain behaves as a visco-elastic 
solid. This visco-elastic mass must be represented in some manner when 
the elastic properties of the booster are calculated. The simplest and most 
straightforward method of accomplishing this is to consider the grain as 
an inert mass, rigidly attached to the case. This method, while it has 
several shortcomings, is in wide use and has been found to yield satisfactory 
results. 

The visco-elastic properties of the grain could be used to provide a more 
comprehensive analysis of the elastic motion. There are several analytic 
models which adequately describe the dynamic behavior of the visco-elastic 
solid (Reference 15 or 16). However, it is generally felt that this area of 
analysis does not need to be considered for study of lateral bending. 

There are several reasons why the visco-elastic properties of the solid 
propellant grain are not used in calculations of the booster elastic properties. 
First, they are found to be relatively unimportant for booster vehicles having 
a reasonable slenderness ratio. The grain structure, in response to stress, 
exhibits a complicated behavior which can be represented as instantaneous 
elasticity, delayed elasticity, and viscous flow. For small stresses occurring 
for short times, the properties could be approximated by considering only 
the range from 500 to 2000 psi at an ambient temperature of 70-80° F, Thus, 
the contribution to the bending stiffness is quite small compared to that 
of the vehicle shell, which is commonly referred to as the solid propellant 
rocket motor case. 

A second consideration is the variable nature of the grain properties them- 
selves, The nature of the approximations whioh oan be used for the model 
to represent the grain would vary depending on the stress level within the 
grain, frequenoy of the application of stress, and temperature, The modulus 
of elasticity is quite temperature -dependent, exhibiting a ohange of roughly 
afaotor of 10 for every 40 degrees F of ohange in grain temperature. This 
property alone makes it oumbersome to desorlbe adequately the solid propellant 
grain motion. This difficulty in analysis, along with the relative unlmportanoe 
of the vlsco-elastio effeots on the mode, has prompted most analysts to omit 
these effeots from the model used to desorlbe the lateral elastio motion of the 
booster. Bending mode tests run by various airframe manufacturers have 
indioated that these omissions do not affeot the adequaoy of the oaloulatlons. The 
above should not be taken to imply that the visoo-elastlo behavior of the solid 
propellant grain is not important in all problems. It does beoome quite Important 
under oertain conditions, particularly in the analysis of the longitudinal modes. 

From this it follows that aspeots Important for the liquid oyllndrioal vehlole 
lateral model are also to be considered for solid boosters. Propellant sloshing, 
of oourse,does not exist, Because of the thioker tank walls of solid boosters, the 
effeots of adapter stiffness and joints are more predominant in the lower frequenoy 
modes and should be carefully examined. 
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3.5 CLUSTERED BOOSTERS. One method for obtaining the higher thrust required 
for large payloads which minimizes some of the penalties in manufacturing costs and 
fuel slosh forces involved in increasing booster diameter is the use of clustered 
boosters. For liquid boosters a peripheral ring of propellant tanks is attached to a 
center tank and the engines are supported on truss members connecting the tanks; 
for solid boosters the motors are attached to a central solid or liquid booster. These 
clustered tank designs destroy axial symmetry and also in some cases planes of 
symmetry. This results in a more complicated lateral model where a number of 
cylindrical tanks are coupled by their elastic connections and in the worst cases must 
be allowed freedom in several directions for an adequate description of vehicle modes. 

For preliminary design it is sufficient to choose approximate planes of symmetry 
and analyze the vehicle for bending modes in pitch and yaw planes using branch beams 
connected to the central core by translational and rotational springs. Simplified 
torsional and longitudinal models will also suffice at this stage. These simple models 
can be used to identify possible problem areas (such as relative modal frequencies) 
and provide design criteria for the connections between tanks. In some cases the 
vehicle characteristics are not, and sometimes even cannot, be known with sufficient 
accuracy to justify any consideration of the effects of non-symmetry. 

A complete analysis (or test) should be undertaken to describe all the primary 
modes of the clustered vehicle. This analysis would provide displacement and 
rotation in two mutually perpendicular planes (preferably principal axis), torsion, 
and longitudinal motion. The model of the tanks for displacement and rotation in 
each of the two planes would be very similar to that discussed for the cylindrical 
booster. Provision must be made to account for the motion of the outer tanks in 
these two directions due to the torsional displacement of the center tank and the 
elastic connections. Longitudinal motions of the outer tanks can couple with the 
bending motion of the center tank; it is also possible that longitudinal motion will 
couple with lateral and torsional displacement. As an example, consider a cluster 
arrangement where the connection at the bottom transmits moment, shear, and 
axial restraints while the connection at the top provides only shear restraint. Then 
it is possible to find a mode where the external tanks are bending, causing moments 
and deflections at the connection to the center tank which will result in longitudinal 
motion of the tank, as well as modes in which the longitudinal motion of the outer 
tanks cause bending of the center tank. The significance of these types of modes can 
only be ascertained from the analysis (or test) and can vary greatly from vehicle to 
vehicle . 

The torsional properties in the model can be represented by the torsional stiffness 
and roll inertia of each tank. The tanks must then be connected by the elastic 
properties of the truss. Formation of the torsional and longitudinal models are given 
in the monographs on those subjects. The complete model for the clustered booster 
then consists of the lateral model in two planes, the torsional model and the 
longitudinal model. These models are then combined through the elasticity and 
geometry of the connections to provide the stiffness and/or mass coupling. 

Analysis (or test) will probably show some planes of symmetry and can also show 
that some of the coupling mechanisms are unimportant for the particular problem to 
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be solved. If this is the case, it is justified, and expedient, to devise a mathematical 
model from the results of the complete anlysis (or test) which represents the modes 
of interest. 

The analytical representation of the clustered booster is more approximate than 
that of the one-dimensional beam lateral model, especially for the condition of coupled 
modes, The tri-axial strain relationships are not completely satisfied, and the 
inertia terms could be poorly represented in the combination lateral-axial motion. 

These effects generally are of second order and as such should not alter the primary 
modes. The elasticity of the connection points probably will require test data for 
accurate values. The numerical techniques employed for the solution of the system 
characteristics have been used for many years in the industry; however, some problems 
in accuracy can be encountered when applied to the clustered booster. The number 
of points describing the system may be compromised for efficient computer operation 
and this vehicle may have modes of nearly equal frequency which will be difficult to 
separate analytically (and experimentally). 

The elasticity of the vehicle is most easily described by a coupled stiffness matrix, 
but solving for the characteristics would then involve inversion of this matrix. This 
inversion may lead to errors because of machine round-off errors. If the inversion 
is successful, round-off errors could still be significant in the iterations for the 
characteristic values, especially for modes with nearly equal eigenvalues. The 
inversion problem can be circumvented by writing the coupled flexibility matrix 
directly. The required transformations can be complicated and lead to errors, but 
with special care this is not insurmountable. Another approach would develop an 
uncoupled flexibility matrix and perform transformations of coordinates to provide 
all necessary coupling in the mass terms. These last two approaches would still 
have possible problems in the iterations on the systems characteristics. 

An approach useful in calculating modes of complex systems is the component 
mode synthesis method. Here the modes of the individual pieces are calculated and 
then the combined modes are obtained from the modes of the component parts. This 
is based on the assumption that significant motions of the individual tanks can be 
described by a small number of modes. If this is true, then the solutions for the 
combined system can be performed in terms of less coordinates. 

Most of the clustered booster work to this date has involved two vehicles — the 
Titan IIIC and the Saturn I (and the second generation Saturn IB) . Titan IIIC is com- 
prised of a center core liquid booster with two attached solid boosters (Figure 2); 
the connections at the bottom transmit axial load, shear, moment, and torque. 

The top connection transmits only shear. Because of the nature of the connections, 
it can be seen that yaw bending and longitudinal coupling can occur and pitch bending 
and torsion is another possible coupling mechanism. Storey in Reference 17, 
develops the coupled flexibility matrices for these two conditions. This method 
encountered difficulty in that the number of stations required for adequate repre- 
sentation of the system with the required transformations exceeded computer capacity. 
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The final Titan IIIC analysis presented in Reference 18 utilizes the mode synthesis 
approach. The longitudinal, torsional, and pitch and yaw bending modes are 
determined for each tank and are then coupled by the elasticity of the connecting 
elements. The influence coefficients for these trusses were obtained experimentally. 

A report giving comparison of analytical and 1/5 scale experimental results is to 
be published. 

The Saturn I vehicle consists of a center LOX tank with eight peripheral tanks 
for alternately LOX and RP-1. These tanks are connected at top and bottom by 
trusses (Figure 3) providing axial, shear, and torsion restraint in both planes at 
the bottom plus moment restraint in the tangential planes. The top connection 
provides similar restraint except for the fuel tanks which do not transmit axial 
load. The trusses are not symmetric with respect to planes of symmetry of 
the tanks, but this effect is small such that planes of symmetry as defined by the 
tanks do not introduce large errors. 

Kiefling ( Reference 6 ) uses a mode synthesis approach for calculation of 
Saturn I modes. Pitch, yaw, and torsion are considered uncoupled and the effect 
of longitudinal propellant vibrations in the outer tanks is coupled with bending. 

A comparison with test data shows very good agreement in frequency, while 
agreement of mode shapes is fairly good. The discrepancy, which is seen in the 
seventh and some higher modes, is in the displacements of the booster center 
tank. Since no control sensors are located in this area, this discrepancy is of 
limited importance for stability and control studies. The mode shape differences 
are due to deflections of the spider beam, the structure connecting the top of the 
booster tanks and the second stage. 

Two higher modes determined analytically were not found during the test. 
Theoretical response in these modes was very small; therefore, the modes would 
be difficult to excite. 

Milner (Reference 19) establishes theoretically the uncoupling of pitch, yaw, 
and torsion modes for a symmetrical clustered booster and investigates the effect 
of minor asymmetry. Results of this study indicate that the effect of such coupling 
on natural frequencies is minor; mode shapes are not presented. 

O'Rourke (Reference 20 ) applies Kron's method to the analysis of an early 
Saturn vehicle. Based on a total of 135 mass points, seven modes of the vehicle were 
found in this analysis. Comparison of these results with dynamic test data is generally 
as good as can be expected, since the analysis was based entirely upon theoretical 
stiffness of the very complex structure. The dynamic data (modes) were obtained 
by a frequency analysis of the vehicle response to a sinusoidal force. It was assumed 
that the vehicle natural modes existed at points of maximum response. 

Bost (Reference 21 ) analyzes the Saturn vehicle using a composite mode method. 
Eight modes are found. Agreement with published test data is only fair. Again, the 
difficulty was of obtaining a structural mathematical model which represents 
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adequately the actual vehicle. 


Republic Aviation currently is engaged in a Saturn I vibration study using Lanczos f 
method of minimized iteration ( Reference 22 ) for solution of the eigenvalue problem. 
The method used for calculating the matrix of influence coefficients is described 
by Meissner (Reference 23 ), while the calculation of the dynamic matrix is shown by 
Berman (Reference 24 ). The complete method is demonstrated on a simple cluster 
beam configuration in Reference 25 , in which seven modes and frequencies of a 

72 degree of freedom system are determined successfully. Preliminary, unpublished 
results from the attempted application to a full scale Saturn including sloshing 
and using 368 degrees of freedom indicate that problems with scaling, round-off 
errors, and excessive machine time exist for the large number of modes desired. 

Whetstone at Lockheed is currently working on an extension of the mode synthesis 
approach to include the effects of spider beam flexibility, but no data has been 
published at the present time. 

Glaser, (Reference 26), constrains the center tank of Saturn I vehicle to oscillate 
in a plane of symmetry; the outer tanks are restricted to motions symmetric to this 
plane. The inner and outer tanks are coupled together through the use of equations of 
compatibility at the connection points. The solution is carried out by a Myklestad - type 
analysis. Lowey (Reference 27) also uses a Myklestad approach on a model whose 
center tank is constrained to move in a single plane; axial motion of the outer tanks is 
permitted. Results from both these analyses are not available at this time. 

Lianis (Reference 28 ) develops a matrix solution of the dynamics problem of a 
four tank booster without center core. The flexibility matrix of the whole unit, with 
appropriate end fixity, is derived. This flexibility matrix together with suitable 
mass matrix is used to derive equations of free vibration in matrix form. The tanks 
are assumed to be similar, but the solution can be modified accordingly for the case 
of non-similar tanks and for other tank configurations. The formulation is general so 
as to furnish any complex mode of vibration. Simple modes, however, can be 
obtained as particular cases of the general problem. 


4/METHODS FOR SOLUTION 


4. 1 FORMATION OF COUPLED EQUATIONS 

The primary purpose of the lateral model is to obtain a representation of the real 
system which can then be represented in mathematical terms. The general approach 
to dynamic problems as given by Equations 1 to 8 requires the formation of the mass, 
stiffness, and dynamic matrices. 

4. 1. 1 STIFFNESS MATRIX. Formation of Equation 8 requires the mass matrix and 
the inverse of the stiffness matrix. Since the inverse of the stiffness matrix is the 
flexibility matrix, one may question a method of analysis beginning with the stiffness 
matrix when a flexibility matrix can be derived directly. The main advantage of the 
stiffness approach is the straight forward manner of deriving a coupled matrix which 
lends itself toward formulation of computer logic capable of assembling a coupled 
stiffness matrix for very complicated systems. Also, with high speed, accurate com- 
puters available, the matrix inversion can usually be accomplished efficiently and 
accurately. Therefore, it is of interest to derive one general computer program which 
will develop the stiffness matrix of various complex systems from simple Input format 
and Instructions. For certain specific problems it may be desirable to develop the 
flexibility matrix. In this section the approach (as presented in Reference 29) for 
developing the stiffness matrix for a lateral model will be given. The flexibility 
approach is presented in Section 4. 1. 2 and a mass coupling technique Is presented 

in Section* 4. 1.3. 

Bending and shear stiffness are identified as points along the structure including 
all mass stations. The nature of the stiffness distributions may Justify stiffness 
definition at an intermediate point between mass concentrations. Thus the model is 
formed from a series of connected massless beams, with mass and inertia concentrations 
located at some or all of the Junctions of the beams. While the stiffness of the vehicle 
may be distributed in a complex fashion, it may be represented with acceptable accuracy 
by a series of straight line segments. This results in giving each beam segment a tra- 
plzoidal stiffness distribution. 

4. 1. 1. 1 Free Element Stiffness Matrix. In the approach to be outlined, the 4 x 4 
stiffness matrix of each element of the beam Is first obtained by inverting the 4x4 
flexibility matrix of this element. (The element stiffness matrices are developed by 
this technique because It is easier mathematically and more accurate. ) The matrix is 
then coupled by constructing a matrix composed of the Individual 4x4 matrices. In this 
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coupled matrix, the terms for common points of adjacent elements are the sum of the 
terms for the individual matrices. With this matrix form, restraints or boundary 
conditions are imposed which will represent the system to be analyzed. The derivation 
of the stiffness matrix is now given. 
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Figure 4. Cantilever Beam in Bending 
The bending stiffness for a value of x is 
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By a similar procedure, the deflection due to shear at x = 0 is 
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Putting (12) and (13) into matrix form 
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where C , C , C , are the coefficients of (12) and(13)> Then, by inversion, 
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The stiffness matrix is then 
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Considering a free element, with the same stiffness as the cantilever beam 
(as shown in Figure 5), the equilibrium equations are: 


P + P 


2 


= 0 


M + P 1 + M = 0 • 


Also from Figure 5, 
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Figure 5. Free Beam in Bending and Shear 
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Equation (15) now becomes : 



Expanding this to include point 2 (by use of the equilibrium equations) : 


(17) 



The stiffness matrix for a free-free element has been obtained. The nature of the 
expressions for Cn, C^2> &nd C 22 such that as EI 2 approaches Eli, expressions 
become indeterminate in form. Each has a limit which is the appropriate coefficient 
for a beam of uniform cross section: 

ii jl 

C 11 3a + c 

c =4 

12 2 a 

C =4- • 

C 22 a 

It has been found that the use of these values when lb, d I < 0. 01 produces 
little or no error, hence, within this range, it is recommended that uniform beam 
coefficients be substituted for the original expressions. 
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4. 1. 1. 2 Axial Load Effects. To determine the effect of axial load on the stiffness 
matrix of a uniform beam element, consider the model in Figure 6. With an axial 
compressive force, H, the moment and shear at any section, x, are given by: 



x 


Figure 6. Beam Element with Axial Load 


M = M 1 + Pj x + H ^ 

V =P - H 

1 dx 



(19) 


d x 


- -(.♦*) 


where 0 Is the slope due to bending and 7 is the slope due to shear. 


Utilizing the relationships t = V/KAG and 


d0 
d x 


M. 

El 


the following is obtained: 


2 

d y _ _ d0 
dx^ d x 



or substituting from (19) 



2 

u y 


M u 
H~ 


„ 2 
P u X 

+ — = + 

H 


2 

u 


y 


l 


( 20 ) 


( 21 ) 
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where 


2 

u = 


H 


(i— M 

V KAG j 


El 


Integrating (21) and applying the equilibrium equations gives the four deformations 
y, , y 0 > and 0o in terms of external forces P , P , M , and M . Expressing 
these relationships In matrix form and Inverting gives the symmetric stiffness matrix 


where 


r > 
M 

1 


K 11 


(sym) 



"l 



p i 


K 21 

K 22 




y l 



< 

M 2 

' 

K 31 

K 32 

K 33 


< 

*2 

> 

(22) 

P 2 

J 


K 41 

K 42 

CO 

K., 

44 


y 2 
L J 




K ll~ r 


sin ul 

■(*-*) 


- I cos ul 


= K, 
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and 


K 21 = y (1 - cos ul) - 


” K 41 " K 43 K 32 


K 22- T *44 = - K 42 


K = iL 

K 31 r 


1 - 


sin ul 


•(>-*) 


r = 2(1- cos 


ul) • “ l ('-gfe ) 31nul • 
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The above matrix is good only for a positive (compressive) load. If the load 
were negative, u would be the square root of a negative number. Therefore, another 
derivation is needed to obtain a stiffness matrix for negative (tension) loads. The 
derivation is the same except that the sign of H is changed and 



is used in place of u . As a result, the elements of the stiffness matrix become 




4. 1.1.3 Local Effects. Locally significant components can generally be included 
in the model by attaching a mass and/or inertia at the appropriate location by means 
of springs representing the elasticity of the mounting structure. The stiffness matrices 
of elements representing translational and rotational springs of Figure 7 can be written 
directly as: 
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Figure 7. Translational and Rotational Springs 




4. 1. 1.4 Coupled Unrestrained Stiffness Matrix. The unrestrained stiffness matrix 
is generated from each beam element, and a 2 x 2 stiffness matrix is generated from 
each spring element. 


Consider a beam element. If an element is attached to a beam, this 
element must either be another beam or a spring, Figure 8. 
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When two beam elements are connected to the same node, coupling with respect to both 
translation and rotation occurs. 

The stiffness matrix layout corresponding to Figure 8 a is given in Figure 9. 

An element in the matrix of Figure 9 gives the magnitude of the force P. or moment Mj 
at node i due to a unit deflection yj or slope at node j. Notice that the two beam 
stiffness matrices overlap at node 2, indicating node 2 feels the effect of both beams. 

An element in this portion of the matrix is obtained by adding the corresponding element 
of the two individual stiffness matrices. 



P = Force 
M “ Moment 

y m Deflection 
8 - Slope 


Figure 9. Stiffness Matrix Layout for Attachment of Two Beams 


The sitffness matrix layout corresponding to Figure 8 b is given in Figure 10. 

The two stiffness matrices again overlap at node 2, but this time at one element. 
Coupling a beam with a spring only affects one degree of freedom, for a spring can 
have but one degree of freedom, (if Instead of a translational spring connected at 
node 2, a rotational spring were connected, then coupling would be at element M 0> 8 0 - ) 
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Figure 10. Stiffness Matrix Layout for Attachment of Beam with a Spring 


4. 1. 1. 5 Reducing the Stiffness Matrix. The stiffness matrix, as developed In the 
preceding may contain coordinates which can be eliminated by consideration of the 
restraints put on a system and the application of boundary conditions to the system. 
Together, the restraints and boundary conditions can be considered to be in two 
classes. These are: 


a. Nodes that have zero generalized displacements but may have non-zero 
generalized forces. 

b. Nodes that have zero generalized forces (no mass or inertia) but may 
have non-zero generalized displacements. 

The equation {P} - [K] M may be rearranged as follows: 



M 



(25) 
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where the P s and M s are forces and moments, respectively, that have the non-zero 
displacement and slopes 0 S and « s . The P 0 and Mo are forces and moments respec- 
tively that have the zero displacements and slopes and « Q . If 



Equation (26) is equivalent to the following two matrix equations 


! F si = M K( 

N - i>'] K f • 


(27) 


(28) 


If there are no boundary conditions associated with the generalized forces 
|f s | , the matrix [a] is the final stiffness matrix [k] and j f s j is a mode shape 
{ \ft) • If there are also boundary conditions calling for zero generalized forces then 
matrix equation (27) is rearranged to become 



( 29 ) 
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where the Ff are the non-zero generalized forces that have the generalized displace- 
ments *f, and the Rf are the zero generalized forces that have the generalized dis- 
placements 

Equation (29) is equivalent to the following two matrix equations 


M Ki + w j» f | . JM 


and 


Therefore, 


[a'J j* f | ♦ [g] j« f ( - )R ( [ = j0( 


and 


I'fi - 


(30) 


K|-[M-[e] [g] -1 [e'j] )f ( [.[K] j* f |. (31) 


In this case (the case of zero generalized forces) the matrix (m 
I s the final stiffness matrix [k] and j^f j Is a mode shape j^j. 


- [•] W [<]) 


4. 1. 1. 6 Elimination of Rigid Body Modes. The foregoing analysis was made for a 
restrained system, one that Is fixed at least In one place. For a structure that Is 
completely free (a free-free system), It can be shown that the eigenvalve problem 
Is one of the form: 


[C] [M] [>] = X[j#| " j*| w Q - jrj (32) 

In the analysis of a free-free system the structure must be fixed temporarily 
at one point. If the structure were not fixed, an external force or moment applied to the 
system would cause the whole system to move uniformly. The solution to the problem 
would then be Impossible. Mathematically this Is represented by a singular stiffness 
matrix [k] . In Equation (32) the term - w q releases the fixed point with respect 
to rranslation, while the term - |t| v 0 releases the fixed point with respect to rotation. 
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The translation and rotation imparted to the fixed point when it is released are w Q and 
v Q . By applying the principles of linear and angular momentum to the system, 
equation (32) can be expressed as a standard eigenvalue problem. 


[ [C] [M] + [TT] + [TR] ] j* j 


X j * j (33) 


where, as explained in the foregoing 


[TT] j * j = » j, j w 0 , 

[TR] j *| * X j, | v 

and 

[ [ c ] [m]+[tt] + [tr]]. 


is (he final dynamic matrix. The derivation of [TTjand [TR] is now given. 
Referring to Equation (32): 


M . = 1. 0 if point i has a translation degree of freedom. 

j! = 0 if point i has a pitching degree of freedom. 

T . - the perpendicular distance from the pitch axis through point zero 
(origin) to point l if point l has a translational degree of freedom. 
(Positive if measured forward of point zero. ) 


r - 1. 0 if point l has a pitch degree of freedom. 
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Applying the conservation of linear momentum 

m o W o + { M 2 | M j * i = 0 (34) 

and the conservation of angular momentum 

l m 0 »o + t T f [“] = 0 < 35 > 

Multiplying equation (32) by JmJ yields 

j » j [m] I* ( - W o |* j [m] j„| - , o j » j' [m] jr j 

* " 2 l»\ M [ c ] [“] | + i • 

Defining 

L A J = j M | [m] [c] [M] 

and substituting from equation (34) 

- mw 0 -w 0 |„j [«][«•] -* 0 jc | [m] Jr j - J LaJ j*(. 

The total mass is 

m t* m <> + 1 <* r [“] )<* ( 
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and the static mass moment about a pitch axis through point zero is 


s= )M|'[M]jr|. 


Therefore 


RL w + S v = 
To o 


- w 2 LaJ | • 


(36) 


Multiplying Equation (32) by | t |[m] yields 

! T ! i M i ! * 1 - w 0 j t M j m | - r o ! T f M I T i 

= * 2 jf j' [M] [C] • 


Defining 


LbJ = j»-([M] [c] [M] 


and substituting from Equation (35) 


"’m " _W I'f MW’' j'f[M]j»j- 

0 0 0 o 

= « L B J . 

The mass moment of inertia of total structure about a pitch axis through point zero is 

m “ T m + 1 T I M )t(. 

T o 
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Therefore 


I 


m 


T 


» + S w 
o o 



(37) 


Solving Equation (37) for v gives 

o 


v 

o 


Sw Q + 


<o 2 LbJ M 



(38) 


Substituting this in Equation (36) and defining 


W = 




S 


2 


gives 


w = 01 

o 




I ~| 

m 

_s 

L B . 

T 1 A 1 

_ w 

j - w L A J 


{+} 


(39) 


Substituting this in Equation (38) gives 


v 

o 


u> 


W 


LaJ 


LbJ 

w 


\*\ 


(40) 


By satisfying the identities 

[tt] )*( * X j,|W 0 
[tr] {+( = X „ o 
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it is seen that the T and R matrices are 



r i n 

[tt] = )m( 

[iLBj - W T LaJ J 


r _ m_ i 

[tr] = |r| 

i 

m 

■ 

< 

i 


4. 1. 2 FLEXIBILITY MATRIX. The Inverted stiffness matrix Is In reality the 
flexibility matrix of the system, i. e. , it expresses displacements in terms of forces: 


|yf = [c] (P| ■ 

If the model is statically determinate, the development of the flexibility matrix directly 
Is much simpler than Inverting the stiffness matrix. However, If the structure is 
indeterminate, the calculation of the elements of the flexibility matrix becomes more 
complex, rapidly becoming Involved and tedious with increasing numbers of redundancies. 
The element, C of the flexibility matrix may be thought of as the deflection at point 
i resulting from a unit load at point j. The principle of reciprocal relations will force 
symmetry of the matrix, reducing the quantity of coefficients to be evaluated. 

The flexibility matrix is formed by, first, developing Individual flexibility 
matrices for each element in the system, considered as cantilevers. 



Note that a) yi and $i are relative to the end considered fixed, and b) Pj and Mi are the 
total loads applied to the end considered free. Consequently, the total deflection value 
is found by transforming yi and from relative coordinates to absolute coordinates: 
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( 43 ) 



The total applied loads, P. and can be considered to be functions of the external 
loads applied to each mass of the structure, expressed by the transformation 



Thus the relationship between total deflection ( and 0j ) and the external loads applied 
to each mass (pj and mi) is developed by substituting the transformed values , 


9 ; 


= [c] 


M, 


M 

H 



M M ft!' 

M [■] ft,!' 

1" C*~| = coupled flexibility matrix . (45) 


For a redundant structure, the influence coefficients are not so readily attained 
and use must be made of an appropriate static analysis such as virtual work on 
Castigliano's theorem (Reference 30). 
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An example of flexibility matrix for a determinant structure is given below: 

Consider the three-mass cantilevered structure shown In Figure 11. The 
quantities and mj are the lateral force and moment applied to the structure at point i. 

If each element is considered to be a cantilever, fixed at the base, the free end deflection 



and rotation relative to the fixed base resultant from the total loads on the free end 
(external loads plus loads transmitted by the preceding element, if any), Pi and Mi, are 
found from the 2x2 flexibility matrix of the separated element: 



These relationships are Independent of the rest of the structure, so the response 
for the entire system can be demonstrated by constructing a flexibility matrix composed 
of the 2x2 matrices of the individual elements: 



The total deflections of the three nodes of the structure, yi and ®i can be written 
in terms of the relative displacements, y^ and J ^ ; 
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3 


y l “ y l + y 2 + y 3 + ‘l *2 + ( ‘l + ‘ 2 ) * 

V *1 + *2 + *3 



In matrix form . 
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Combining Equation (46) and (47) yields, 


1 0 1 l i 1 (4 + , 2 ) »*i 

0 10 10 1 »6 l 9$ 1 


0 0 10 1 1 

mt 

0 0 0 1 0 1 

0 0 0 0 1 0 

0 0 0 0 0 1 


«2 i9 2 


*h "2 


“1 

M 1 

(‘V'i 

"0 ("l* 



» *2 


0 

0 

i<2 

i 0 


•«3 

I U* + l *U*\ 



The total load at any element and can be related to the applied loads at each point by 



or, In matrix form 



( 49 ) 
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which is equal to 
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4. 1. 3 TRANSFORMED MASS MATRK. Another alternative technique for forming the 
dynamic matrix Is to transform the coordinate system from the absolute to the relative 
sense. The equations of motion formed previously consider the displacements of the 
respective coordinates to be referenced to a fixed point, or neutral position. The 
displacements may also be expressed 1 relatively - referenced to an adjacent coordinate. 
Furthermore, the Inherent relationship between the displacements In absolute terms, 
y , and the displacements In relative terms, y , Is readily expressed by a 
simple transformation matrix: 


jyj = [t] jy J 


For example, the simple cantilevered beam In Figure 12 has six coordinates. The 
relative -absolute relationship Is given by: 


y = y + y + y +1 0 + (l + 1^0 

*1 y l y 2 y 3 12 VI 2) 3 


v— y+V+l 0„ 

y 2~ y 2 2 3 1 


,y 3 = y 3 , 


j = 6 +0 +0 , 

1 1 2 3 ’ 


1 = 0+0 
2 2 3 


} = 0 

3 3 
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3 


H 


'2 ^ \ *1 


o 


«2 


-o : O 


Y 2 y x 


Figure 12. Three Mass Cantilevered Beam 


Thus, 


r 

y i 


1 

1 

1 

0 

>2 

(■**'*) 


C* 

y l 

y 2 


0 

1 

1 

0 

0 

*3 


y 2 

y 3 


0 

0 

1 

0 

0 

0 


y 3 

< 

fll 

> = 

0 

0 

0 

1 

1 

1 

< 

> 

• 1 

H 


0 

0 

0 

0 

1 

1 


h 

Jl 

: 

0 

0 

0 

0 

0 

1 


J 3 J 


|y|= [T] jy|. 


(51) 


The kinetic and potential energies of the system can then be written (In matrix notation) 
in terms of relative coordinates as. 


2KE -lyi'[T]'[M] [ T ]|*| 

2PE - |y T [c" 1 ] | yj. 
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where C is the uncoupled flexibility matrix consisting of 2 x 2 matrices for cantilevered 
elements. Using La Grange's equation 


d / 3 ke \ + aPE + aw _ 

dt l dy i J dy { ay L 


the equations of motion become 


[t] [m] [t] j’y | + [ c" 1 ] jy[- 0 


and the dynamic matrix is 


[c][ t ][m][ t ] ])♦!- xj+| 


(52) 


(53) 


(54) 


This approach is very similar to the method in the previous section with the 
transformation of coordinates coupling the mass matrix Instead of the flexibility 
matrix. Note that the modes, j^j , are in terms of relative coordinates and 
must be premultiplied by the transform matrix, [t] , to obtain absolute vectors. 

4. 2 SOLUTIONS FOR CHARACTERISTICS 

Formulation of the equations of motion of dynamic systems results in a linear 
differential equation for the continuous exact solution or a series of differential 
equations for the approximate solutions. For lateral vibration usually only the lower 
modes are of any significance and therefore the approximate solutions are of practical 
importance. Two methods are used to describe the system in these approximate 
solutions: (1) the system is divided Into a finite number of segments connected by 
massless stiffness and (2) the system Is described in terms of assumed functions. 
Solving for the characteristics of the resulting equations can be categorized into 
three groups. These are: (1) energy methods, (2) solving the differential equation, 
and (3) solving the Integral equation. The equations are in the general matrix form: 

-w 2 [mJ j z j + [k] jzj = 0 (differential equation) (55) 

or 

jzj= jz | .. (Integral equation) (56) 
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The most general solution of Equation (55) Involves expansion of the determinant and 
solving the polynomial equation. This procedure is adequate for simple systems and 
up to four degrees of freedom can be solved easily. Methods for higher order systems 
have been developed in References 31 and 32. However, since only the lower modes 
are important some approximate methods have been developed which obtain these modes 
and frequencies with sufficient accuracy. 

4. 2. 1 MATRIX ITERATION (STODOLA AND VIANELLO METHOD). The matrix 
iteration technique is essentially the Stodola and Vinanello method in matrix form. 

The integral equation Is 


| z |= ” 2 [K]- 1 [m] M ’ (57) 

It can be seen that « £ M ] j Z n _^ | , is the load associated with a assumed mode 
shape, | , vibrating at a frequency w . The deflection resulting from this 

load or the next approximation of the mode shape is obtained by premultiplication by 
the influence coefficients, or the inverse of the stiffness matrix. Since w 2 [ S a constant 
it can be assumed to be unity and the equation becomes 

j z n i = [k ]" 1 M I Vi I ' < 58 > 


where j Z*j is a vector normalized on its largest element. Successive iterations of 
Equation (58) continue until jz n *| has converged, i.e. , every element in jz n *j 
satisfies 



The frequency of the mode is then 


w 


2 



(59) 
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where Z n j Is the largest element of the unnormalized vector jz n j . Proof of convergence 
is given in Reference 33. 

For most practical problems it is necessary to obtain more than the first normal 
mode. This is accomplished by applying the condition of orthogonality as a means of 
purifying the assumed higher modes of lower mode components. The orthogonality 
condition requires that 

Ki' [“] I M- °- (60) 


Therefore, modification of the dynamic matrix to satisfy this condition will allow 
extraction of the next mode. This is expressed in matrix form as: 


[D 2 ]-[dJ [[I] -[a]] (61) 

where [i] is a unity matrix and [ A ] is a matrix of zero’s except for the row (or 
rows) imposing the orthogonality. This row is composed of the elements Z. m* 
normalized on its largest element. This element locates the characteristic row and 
the row matrix is a good approximation of the characteristic row. As an example , 


[z' mj - Lz n ■»! Z 21 “ 2 Z 31 ” 3 Z 41 m 4 J 


If Z m is the largest element, then 
31 3 


[A] 


0 

0 


Z 11 m i 
Z 31 m 3 


0 

0 


Z 21 m 2 
Z 31 m 3 


0 

0 


0 

0 


Z 41 m 4 
Z 31 m 3 
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Once £ A J is obtained then D 2 can be formed and iterations proceed to determine the 
next mode and frequency. 

For improved accuracy and subsequent deflation of the D matrix it is recommended 
that iterations be performed on the row matrix, [X^M J , by post multiplication by the 
dynamic matrix. This removes errors in roundoff in the preceeding fundamental mode 
and provides a more accurate characteristic. Thus, 


L z j M J [°] = L z - M J 


where [ ] Indicates normalization on its largest element. L Z? M J is then normalized 

on its largest element and placed in the [ A J matrix. If the characteristic row is found 
by this method, the deflated dynamic matrix for successive modes can then be obtained 
by 


N- [° 2 ] [w - [ a 2 ]]' 

If the deflated matrix is obtained only from the conditions of orthogonality, then 

[°n] ’ [ D l] [M - [AJ] . (62) 

where[ A n ] satisfies orthogonality between first and n, second and n, , and n-1 

and n modes. 


An alternate method for deflating the matrix is given in Reference 34 and is 
particularly useful in application with automatic computers. This deflation method 
gives 

[ d 2 ] ‘ [ D i] ■ l z i! N [M] < 63 > 

“1 
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where {Zj} Is the column of elements of the first mode normalized so that 
[ZjJ [M] {Z^} ■ 1. For each succeeding mode the matrix Is modified by subtracting 
the triple matrix product for the proceeding mode from the proceeding modified dynamic 
matrix 


[*>»] 


P 


n-l_ 


JL 

2 

i) 

n-1 


« 


{Vi} L z -iJ [M] 


( 64 ) 


4. 2. 2 HOLZER-MYKLESTAD METHOD. Tne Holzer-Myklestad mechod Is extremely 
suitable for obtaining solution ox the differential equations of beams represented as con- 
centrated masses connected by massless stiffness elements. A principle advantage of this 
method is that eacn mode is obtained Independently and therefore all modes are as accurate 
as represented by the system analyzed. 

The technique proposed by Holzer was originally developed for torsional systems. 

It is equally applicable to any close-coupled system, i. e. represented as a spring 
mass system. A frequency is assumed and one element of the mode shape vector is 
taken as unity. Through equilibrium equations, the other elements of the vector are 
calculated. There will be one additional equilibrium equation left after the last element 
of the mode shape vector Is evaluated, which will be satisfied only if the assumed 
frequency is correct. Successive frequency trials are made until the proper value 
is obtained. The process is facilitated by plotting net unbalanced force on the final 
mass against frequency. At the natural frequencies of the system, the curve of net 
force passes through zero. 

As an example, consider the simple system of two masses and two springs shown 
in Figure 13. The governing equations are 


m. X. + 

K 

X. - 

K, X n - 0 

1 1 

1 

1 

1 2 

m 2*2 - 

K 1 

X + 

(V K 2) X 2 -° 
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Figure 13. Two Degree of Freedom Spring-Mass System 
The assumption of harmonic motion reduces these equations to 



(60) 


(67) 


The magnitude of the left hand side of this equation can be plotted against frequency. 
It may be facilitated by calculations made In tabular form: 










When the residual Is zero, the correct frequency has been obtained. 


Branched systems will necessitate temporarily expressing some elements of the 
mode shape vector in terms of others, but in all systems, numerical values for the 
entire mode shape will be determined upon solution of the next to last equilibrium 
equation, leaving the last equation for frequency evaluation. 


As an example consider the simple three mass system shown in Figure 14. 
-A * (K, + K 3 )X 1 -K 1 X 2 -K 3 X 3 = 0 

”2*2 ' K 1 X 1 + ( K 1 + K 2 * K 4 ) X 2- K 4 X 3'° 

m 3 X 3 - K 3 X 1 - K 4 X 2 + ( K 1 + V X 3“ ° 

X, 



x„ 


Figure 14. Three Degree of Freedom Spring-mass System 


Introducing harmonic motion results In 

( K 1 + V " I "l) X l- K l X 2 - K 3 X 3 

-KjXj + (Kj +K 2 + K 4 - « 2 m 2 )X 2 - K 4 X 3 = 0 
- K 3 X 1 - K 4 X 2 + ( K 1 + K 4-« 2 m 3 ) X 3 '» 
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Let Xj = 1, and from the first equation 


V 


K 1 + V " m i~ K l X 2 

K 


From the third equation 


- K 3- K 4 X 2 


+ ( K 1 + K 4 - “ »») 


K 1 * K 3 - “ m i~ K 1 X 2 

K„ 


or 


K 3 + ( Kl + K 4 - w m 3 


\ + K 3 “ “ 


m. 


K_ 


K. . 

K 4 + F 3 ( K 1 + K 4- " X) 


x 2 =o 


X 2 = 


_f K l + V ^ 3 )K * K 3 - “ ‘ K ; 

K 2 K 4 +K 1 (\ + K 4- « m 3 ) 


2 

3 _ 


= A 


(68) 


Hence 


+ Kg - (•) m^ - A 


= B 


( 69 ) 


And finally, from the second equation 


V ( K 1 * K 2 + K 4- " m 2 ) A - K 3 B=0 - 


( 70 ) 


The three frequencies which satisfy this equation are the natural frequencies awl are 
determined by tabular methods as In the previous example. 
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In reality, the Mykiestad method is the Hoizer method extended to include far- 
coupled systems, l. e. beam systems. The slope and deflection at one end of a beam 
segment are expressible in terms of the slope and deflection at the other end, the 
loads at the other end and the flexibility of the beam segment. The boundary conditions 
at one end of the beam give an initial set of values. Two of the four conditions will 
be unknown. One is set to unity and the interior values are determined in terms of 
the remaining unknown. 



Figure 15. Free-Body Diagram of Vibrating Beam Segment 



The quantities (64),(40),(0j),(00) are the flexibility Influence coefficients 
of the beam. If the right end of the beam segment Is the end of the total beam (let n = 1), 
two of the four valables at that end can be determined from the end conditions: for a free 
end. Mi is zero and Vi is equal to the inertia force mi w 2 y . Of the two remaining 
conditions, one Is assigned a value of unity and analysis Is carried out in terms of the 
other. For the free end, Xi is assigned a unity value (making Vi ■ mi w2). The four 
variables at point two may be evaluated as: 
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y 2 = 1 " m i “ 2 ( a< ) • l x 9 i ■ 
*2 =*1 " m i U (**) ’ 


,, 2 . 

M 2 = m l y 1 , 

2 * 
V 2 = m i w " m 2 


1 - 



If an arbitrary value fory2 Is assigned, each of the four variables at point two 
reduces to an expression Involving a constant and a coefficient times 8j_. 

The variables at the next point on the beam may be expressed in terms of those 
just obtained by using Equation (71) and will also be in linear terms of The 
evaluations may thus be propagated across the beam, with the ultimate result that 
the four variables at the far end of the beam are expressed as linear functions of 

From the end conditions at this point, the values of two of these variables will 
be known. One of these known quantities Is used to evaluate which is then 
substituted into the expression for the second known variable. For the correct value 
ofw 2 , the expression will produce the known value; otherwise, an error term or 
residual will be obtained. Several trial values oij^ are made, repeating the pro- 
pagation process each time. If the resultant residuals are plotted agalnsty^, the 
extraction of the proper values ofy^ [ s facilitated by extrapolation of the curve. The 
natural frequencies occur where the curve crosses the axis. Higher frequencies 
are obtained independently, thus are Independent of any inaccuracies which may 
exist in calculated lower modes. 


4. 2. 3 ENERGY METHODS. Lord Rayleigh's method of evaluating the fundamental 
frequency of a system is based on the principle of conservation of energy. At the 
maximum deformation of the system vibration in its fundamental frequency, all the 
energy of the system is in a potential energy from ^PE = J~El ^ d2y ^ 2 dx \ . 

But at the instant the system passed through the equilibrium position, its energy is 
entirely in kinetic form ^KE = ^/my 2 dxj . If energy is conserved, the maxlmums 
of those two values may be equated. A deflection shape Is assumed and for harmonic 
motion, y = U sin wt. , and hence, 


y = « U cos wt 


2 J" EI f ) dx “ \ | m U 2 y 2 dx , 


(72) 


and thus 


» 2 = / ei (^-) 2 d: 


(73) 


/ 


m U ^ dx 
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It can be shown that close approximations to the fundamental frequency are 
obtained even if the assumed deflection shape is not very close to the true shape. 

The Rayleigh method was generalized by Ritz to give more'accurate 
values for the frequencies as well as to give estimates for several mode 
frequencies at one time. Basically, Ritz suggested that the assumed deflection 
curve be expressed as the sum of several functions in the form 


y = 2 f (x) Q n . 

n a 


(74) 


The more functions and constants Introduced the more accurate will be the value of the 
fundamental frequency. Also if n functions are introduced estimates on the first n mode 
frequencies will be obtained. Having expressed the assumed deflection curve in terms 
of n functions f (x), the kinetic and potential energies (omitting shear & rotary inertias) 


are expressed as 


2 KE =Jra (x)y (x, t) dx 

(75) 

2 PE =/EI (x) dx . 

(76) 


Substituting Equation (74) into Equations (75) and (75), 


n n 



2 KE = 2 2 

i=l j=l 

m u q i q j 

(77) 


n n 

2 PE = 2 2 

i=l j=l 

K u "l “i 

(78) 

where 

K n =/ EI 77 [ 

V x >] i? [*,«] dx 


and 

m i)=/ m 772 [ 

V x) ] S [ f j < x >] dx ■ 
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Substituting Equation (77) and (78) into La Grange’s Equation gives the differential 
equations 


& V£Vi-°- 


(79) 


Putting as a solution q^ = 7} sin u»t 
/ 2 

( K ij “ w m I} ) = 0 (80) 

or a set of algebraic equations equal In number to the number of unknown coefficients. 
The determinantal equation for the set will yield frequencies for the first n modes of 
vibration and eigenvectors which when multiplied by the assumed functions give 
approximate mode shapes. 

The result gives good accuracy In frequency but poor agreement with 
mode shape. Using the mode shape obtained as a second approximation 
improves frequency accuracy and greatly improves mode shape. 

Also, expressing the strain energy in terms of inertial loading rather 
than the assumed deflections increases the flexibility of this method in that 
shear deformation and rotary inertias are then easily included. 

4. 2. 4 MODAL QUANTITIES. Solutions to the characteristic equations give 

the reciprocals of the squares of the circular frequencies and also the mode 
shapes of the restrained system. The linear frequencies of vibration are 
obtained from the circular frequencies. If each mode shape, {^} n is considered 
to be the nth column of a matrix [♦] of all the mode shapes then, [♦J 4 * * 7 [M] [♦] 

is an orthogonality check of the mode shapes. The diagonal element ^ [tf*] 7 [M] [$] j 

is called the generalized mass of the system for the mode n. (The generalized 
mass of the system may be considered to be a measure of the kinetic energy of 
the system). In this discussion the vector was stated to be mode shape of 
the structural system. For a restrained system, is the complete mode 
shape. For a free-free system the displacements due to w Q and y 0 must be 
added to obtain the complete mode shape. Also, for a free-free system, the 
generalized mass must be modified to Include the contribution of the temporarily 
fixed point. 
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The slope of a node of a beam as computed by the foregoing Is not the total slope 
of the beam at the node. The slope computed Is the slope due to beam bending. To 
obtain the total slope, the shearing slope of the beam must be added to the slope due to 
beam bending. This is done by 



w- 

•} + H 

(81) 

where 

M = 

[ss] {fb} 


and, 





| 0 j = total slopes of a beam element , 

| g | = bending slopes of a beam element, 

| y | = shearing slopes of a beam element, 

{fb[ = internal forces on a beam element , and 
{ss} = shear - slope matrix . 

The matrix [ss J Is a diagonal matrix whose elements are l/KAG at the nodes of 
a beam element. 

The generalized forces of the restrained system may be obtained from either 
Equation (31) or Equation (27) depending upon whether or not there are any zero 
generalized force boundary conditions. Rather than use Equation (31), another 
equivalent relationship is used. Extending Equation (31) 

{ F f } = [K] [ * ] = W 2 [M] j* | . (82) 


Equation (82) is used instead of Equation (31) when calculating the generalized forces 
of the restrained system. Using Equation (30) the generalized displacements that have 
zero generalized forces are calculated. For a free-free system Equation (30) is modified 
slightly. In its stated form Equation (30) would compute generalized displacements with 
respect to the temporarily fixed point. For a free-free system the analogy to Equation 
(30) is : 
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{'«} -[ g ] *[ e ] {" w o + + {"Vo + {'}«o (83) 

/ / 

where p, T , w 0 and v Q have previously been defined, and n and t have the same meaning 
as m and r. However, and r refer to the mode shape, , and j/and r' refer 

to the generalized displacements that have zero generalized forces, | Pf } . The 
generalized forces that have zero generalized displacements are obtained from Equation 
(28). The internal forces and/or moments (shears and/or bending moments) of each 
element are obtained from : 

(p}= [*]{*}• 


4. 3 MODE SYNTHESIS ANALYSIS 

The complicating aspects of the clustered booster vehicle have, as a net effect, 
the requirement of large numbers of coordinates in the model. The resultant size 
of the governing equations may well be so large as to overwhelm the best of analysts 
or computers. The technique of modal synthesis Is a process whereby the dynamic 
characteristics of the several components of the system are calculated separately, 
and then brought together to evaluate the dynamic characteristics of the entire 
system. This technique is also applicable to the analysis of staged vehicles in that 
the natural vibration characteristics of the Individual stages can be computed sepa- 
rately and combined via modal synthesis to determine the characteristics of the entire 
vehicle. 

The entire concept of modal synthesis is founded on the principal that a component 
is completely (or adequately) represented by its primary modes and if this is true, then 
the connections of these components at their Interfaces can be described in terms of the 
component modal quantities. A set of equations can be written with coordinates in terms 
of component mode shapes and amplitude factors. Solving these equations for the amplitude 
factors and multiplying by the component modes will give system modes. 


The division of the system into components can be done by removal of connecting 
flexibility elements. The parts thus formed will constitute free system l.e. , not attached 
to a fixed point. It is assumed that the components will be defined in such a way that 
determination of their individual dynamic characteristics does not become a complicated 
task and consequently, the techniques previously discussed regarding the analysis of 
simple systems may be utilized to provide these dynamic characteristics. 
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For the free components, complete description of motion requires the rigid body 
translation and rotation modes and the elastic modes. Thus, If the component Is In beam 
form, the lateral deflection of point Is given by: 


y l =Y + L 1 ® + Z n *ln "n 


where 

y^ Is the lateral displacement of point l , 

Y Is the rigid body lateral translation of the component c. g. , 
Lj Is the distance between the component c. g. and point l , 

0 Is the rigid body rotation of the component c. g. , 

71 n Is the appropriate weighting factor applied to mode n , and 
Is the lateral displacement of point l in mode n , 

The rotation of point l Is given by 


1 = 0+1 
l n 


V 

n 


where 

$ Is the angular displacement of point l , and 
Is the slope at point l In node n . 

A simple component beam will have Its coordinate motions expressed In matrix 
form as : 
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More generally 


|y| = [Ts]jv|. 

For a system composed of several component beams 


r > 

y w 

y< 2) 


TS^ 1 ) 

T g( 2 ) 


r 

, <2 » 

< ; 

• (S) 
y 

I! 

9 

• j 

9 

9 

• 

* 

_ '''■•TS (S > 

< 

; > 
V) 

J 


The kinetic energy of the system is then 

* f / , 

2 KE = t| (TS) M (TS) H . 


and 

jsL ( *ke\ , 
dt - ( TS ) M (TS) n 

where M is a diagonal of system masses. 

Note that (TS) M (TS) is a diagonal of the generalized masses of modes used in the 
representation 

n (!) \ 

? (2) ( ( 86 ) 


(TS) M (TS) n = 


m 


(1) 


m 


( 2 ) 


m 


(S) 
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Similarly, the potential energy In the component modes can be expressed as 


2 pe = n '(TS)'k (TS) n 


and the energy In the connecting elements can be expressed as I, K „A where K ck 
Is the stiffness of the connecting element and'A^ Is the deformation of that element. 

If y r and 0 r describe the coordinate displacements of the components, then 


{A} - 



(RC) (TS) n 


(87) 


where (RC) Is a transformation matrix expressing the relation between and the 
component displacements. The potential energy la given by 


2 PE - n' (TS)' (RC)' K (RC) (TS) 1 

0 

and the total potential energy of the system 


2 PE - ii 
dPE 


(TS)' (RC)' K (RC)' (TS) n + n 

C 

[(TS)' (RC)' K c (RC) (TS) + 


(TS)'K (TS) i) 
(TS)' K (TS)'] n . 


It is shown In Reference 35 that 


(TS)' K (TS) = > 


Therefore, the complete set of equations Is 


fm 1 {M}+[h- 2 mJ + [(TS)' (RC)'K c (RC) (TS)]] (n) * 0 (88) 
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There are other ways to develop a modal synthesis technique for complex systems 
which are quite valid. But the approach presented above enables the analyst to evaluate 
a minimum of physical parameters, then manipulate them with matrix techniques, which 
are readily adaptable to computer programming, Into the final form of the governing 
equation: the parameters to be obtained are just the component modal properties (mode 
shapes, generalized masses, frequencies), the coordinate displacements In terms of the 
rigid body displacements and modal weighting factors (the [TS] transformation matrix) 
and the deformations of the connecting elements In terms of the coordinate displacements 
(the [RC J transformation matrix). The final governing equation may be formed Immediately 
from these data. Solutions for these equations can be obtained by the methods given In 
Section 4. 2. 

These roots will be the frequencies of the coupled system, and the vectors \ n 1 
when premultlplled by [TSjwill give the system modes. 


To demonstrate the modal synthesis technique, consider the spring-mass system 
of Figure 16. 



Figure 16. Three Degree of Freedom Spring Mass System 


The close coupled system Is used for simplicity. By removing spring element K , 
the system Is divided into two smaller components. It can be shown by elementary 
dynamic analysis that the dynamic characteristics of the two component parts are: 
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Component 1 Component 2 
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The displacements are expressed In terms of the component modes as: 


or 


x i 

= Xi * 

*11 

’l 

X 2 

= x i * 

*21 

n l 
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Ths deformation of the connecting element Is given by 



which can be simplified to 



The contribution of the connecting element to the stiffness matrix Is 
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or 



The contribution of the component Is the diagonal matrix 
0 _ 2 


K. 


m i +m 2 


m 


2 J 


K 


These two are added to form the final stiffness matrix. 

The mass matrix Is formed directly from the component mode generalized 
masses, 


m i + m 2 


m. 


m. 





The governing equation Is thus: 


72 


♦ 



To examine the accuracy of this technique applied to this particular model, 
numerical values may be applied to the system parameters and solutions obtained from 
the above equations and also from equations derived In the classical manner. 

Let M 1 * 10 slugs (143 x 106 dynes), M 2 - 20 slugs (286 x 106 dynes) Mg - 30 slugs 
(429 x 106 dynes), K x - 15000 lb/ft (224. 8 kg/cm), K 2 - 12000 lb/ft (179 kg/cm), K 3 - 
10000 lb/ft (149 kg/cm). The eigenvalues and eigenvectors resultant from the two sets 
of equations are as follows: 


Eigenvalue 1 
Eigenvector 1 
Eigenvalue 2 
Eigenvector 1 
Eigenvalue 3 


Classical Equation 

7.6651127 x 10" 3 

1.00000 } 
0.91303 > 
0.60578 ) 

1.1096561 xlO -3 

1.00000 } 
0.39921 V 
-0. 95138 ) 

3. 9189681 x 10 _4 


Mode Synthesis Equa. 

7.6651123 x 10" 3 

1.00000 ) 
0.91303 [ 
0.60578 ) 

1.1096558 x 10“ 3 

1. 00000 ) 
0.30921 1 
-0. 05138 ) 

3.9189677 x 10" 4 



( 1.00000 } 

( 1. 00000 } 

Eigenvector 3 

< -0.70113 > 

<-0.70113 1 


* 0.15423 ) 

* 0. 15423 j 
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